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O I In nucleosynthesis three possible paths are known to bridge the mass gaps at A = 5 and 

(^ ' A = 8. The primary path producing the bulk of the carbon in our Universe proceeds via 
O ■ the triple-alpha process ^He(2Q;,7)^^C. This process takes place in helium-burning of red 
giant stars. We show that outside a narrow window of about 0.5 % of the strength or range 
of the strong force, the stellar production of carbon or oxygen through the triple-alpha 
process is reduced by factors of 30 to 1000. Outside this small window the creation of 
^ I carbon or oxygen and therefore also carbon-based life in the universe is strongly disfavored. 
C ■ The anthropically allowed strengths of the strong force also give severe constraints for the 
sum of the light quark masses as well as the Higgs vacuum expectation value and mass 
K^ , parameter at the 1 % level. 

^ 1. INTRODUCTION 

Few-body methods are used in nuclear astrophysics for the determination of thermonu- 
clear cross sections and reaction rates predominantly for nuclei with mass numbers up to 
about A = 12. These reactions are of relevance for primordial nucleosynthesis, i.e., the 
production of nuclei up to A = 7 in standard Big-Bang model. In the inhomogeneous Big- 
Bang nucleosynthesis model nuclei with A > 7 can also be produced. Few-body methods 
play also a role in the determination of reaction rates for the pp-chains in main-sequence 
stars up to A = 7 and helium burning from A = 4 to A = 12. Finally, in the alpha-rich 
freeze-out occurring in supernovae, reaction rates for nuclei from A = 4 to A = 12 are 
studied using few-body methods. 

We want to discuss the bridging of the mass gaps at A = 5 and A = 8 in nucleosynthesis 
through three-body reactions (Fig. |T|). The following paths via three-body reactions can 
bridge these mass gaps: 
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Figure 1. Possible reaction paths bridging the mass numbers A = 5 and A = 8. 



1. Reaction path ^He(2n,7)^He(a,n)^Be(a,n)^^C: Under representative conditions the 
reaction path ^He(2n,7)^He(a,n)^Be(a,n)^^C can be neglected in known astrophys- 
ical scenarios, because the fragile ^He is destroyed very effectively by photodis- 
sociation ^He(7,2n)^He ||l]-0|. The cross sections and reaction rates for the reac- 
tion ^He(2n,7)^He and the reverse photodissociation have been calculated using an 
a+n+n approach 0. 

2. Reaction path ^He(an,7)^Be(a,n)^^C: This reaction path is dominant in the so- 
called alpha-rich freeze-out occurring in type II supernovae [jl],^-0. The cross sec- 
tions and the reaction rates for ^He(an,7)^Be and the reverse photodissociation have 
been calculated in a semimicroscopic model and compared with the experimental 
data [0]. 

3. Triple-alpha process ^He(2Q;,7)^^C: The bulk of the carbon in the universe is pro- 
duced through this process in the helium-burning phase of red giant stars. The 
reaction rate for the triple-alpha process has been determined recently in a micro- 
scopic 12-nucleon model |P^IO|. 

In the following section we will present the calculation of the reaction rate for the 
triple-alpha process in the microscopic 12-nucleon model ^He(2a,7)^^C. Sect. 3 will be 
devoted to variations of the nucleon-nucleon (N-N) force and its consequences for the 
stellar production of carbon and oxygen in the universe. In Sect. 4 we discuss the con- 
straints for anthropically allowed parameters of elementary particle physics like the light 
quark masses or the Higgs vacuum expectation value and mass parameter. In the last 
section a short summary is given. 



2. REACTION RATE FOR THE TRIPLE-ALPHA PROCESS 

The formation of ^^C in hydrogen burning is blocked by the absence of stable elements 
at mass numbers A = 5 and A = 8. Opik [|n|] and Salpeter |]T^ pointed out that the 
lifetime of '^Be is long enough, so that the a + a ^ ^Be reaction can produce macroscopic 
amounts of equilibrium ®Be in stars. Then, the unstable ^Be could capture an additional 
a particle to produce stable ^^C However, this so-called triple-alpha reaction has very 
low rate since the density of ^Be in the stellar plasma is very low, because of its short 
lifetime of 10~^^ s. 



Hoyle (cited in [|T^]) argued that the triple- alpha reaction cannot produce enough carbon 
in a non-resonant way in order to explain the measured abundance in the Universe, 
therefore it must proceed through a hypothetical resonance of ^^C, thus strongly enhancing 
the cross section. Hoyle suggested that this resonance is a J'^ = O"*" state at about 
6 = 0.4 MeV (throughout this paper e denotes resonance energy in the center-of-mass 
frame relative to the three-alpha threshold, while F denotes the full width). Subsequent 
experiments indeed found a O"*" resonance in ^^C in the predicted energy region by studies 
of the reaction ^''N(d,a)^2C ^!^ and the /^'-decay of ^^B |5|. It is the second 0+ state 



(0^) in ^^C. Its modern parameters, e = 0.3796 MeV and F = 8.5 x 10 ^MeV 0, agree 
well with the old theoretical prediction. 

2.1. Microscopic 12-body model 

The astrophysical models that determine the amount of carbon and oxygen produced 
in red giant stars need some nuclear properties of ^^C as input parameters. Namely, 
the position and width of the 02" resonance, which almost solely determines the triple- 
alpha reaction rate, and the radiative decay width for the 0^ — >■ 2^ transition in ^^C We 
calculated these quantities in a microscopic 12-body model [p|-[lO[ . 



In the microscopic cluster model it is assumed that the wave functions of certain nuclei, 
like ^^C, contain, with large weight, components which describe the given nucleus as a 
compound of 2-3 clusters. By assuming rather simple structures for the cluster internal 
states, the relative motions between the clusters, which are the most important degrees 
of freedom, can be treated rigorously. The strong binding of the free alpha particle ^He 
makes it natural that the low-lying states of ^^C have 3Q;-structures [|T^. Therefore, our 
cluster- model wave function for ^^C looks like 

^"^ = E ^{$"$"$%tli(pi' P2)}- (1) 

ll,l2 

Here A is the intercluster antisymmetrizer, the $° cluster internal states are transla- 
tionally invariant Os harmonic-oscillator shell-model states with zero total spin, the p 
vectors are the intercluster Jacobi coordinates, li and I2 are the angular momenta of the 
two relative motions, L is the total orbital angular momentum and [. . .] denotes angular 
momentum coupling. The total spin and parity of ^^C are J = L and tt = (— 1)'^+'^, 
respectively. 

We performed calculations using the Minnesota (MN) and modified Hasegawa-Nagagta 
(MHN) N-N forces. These two forces give the best simultaneous description of the ^Be 



ground state and the O2 state of C |10], in agreement with the experience gained in the 



cluster-model description of the structure and reactions of various light nuclei [|1^]. The 



MN-force is given by [|T9|,|0 



V^At) = ( V^{r) + ^(1 + Pr^V2{r) + ^(1 - P^^W^ir) ) ( 1^ + 1(2 - u)P[^ )+yCoui.(^)^(2) 



where P^ and P'^ are the space- and spin-exchange operators, respectively, u is the ex- 
change mixture parameter, r = \rj — ri\, and V^°^^' is the Coulomb force between the two 



nucleons. The MHN force is given by [^ 



V^Ar) = E (W^fc + MkP[, + B,P^^ + H,P[^Vk{r) + V, 



Coul. 



(3) 



fc=i 



The parameter Wk is the Wigner parameter, P^^-, P^^, and Pj'^- are the space-, spin-, and 
isospin-exchange operators (Majorana-, Bartlett-, and Heisenberg-operators) , and M^, 
-Bfe, and H^ are the corresponding exchange mixture parameters. 
The spatial parts of the MN- and MHN-force have Gaussian forms 



Vfc(r) = Vofcexp 



\rokJ 



k = 1,2,3, 



(4) 



where Vqa; and r^k are the strength and range parameters of the potentials, respectively. 

2.2. Triple-alpha reaction rate 

The resonant reaction rate for the triple-alpha process proceeding via the ground state 



of ^Be and the 0^ resonance in ^^C is given approximately by [p2 



rso 



3-^N' 






exp 



knT 



(5) 



where Ma and Na are the mass and the number density of the alpha particle, while h and 
k are Planck's and Boltzmann's constant, respectively. The temperature of the stellar 
plasma is given by T. 

We calculated the resonance energy e of the 02" state in ^^C, relative to the 3a-threshold, 
and the 02" -^ 2^ radiative (E2) width F^. The 02" state is situated above the 3a-threshold, 
therefore for a rigorous description one has to use an approach which can describe three- 



body resonances correctly. We choose the complex scaling method [^ that has already 
been used in a variety of other nuclear physics problems, see e.g. |T7| , p4| , p5| . 

As a first step of our calculations, we fine tune each N-N force (by slightly changing their 
exchange-mixture parameters) to fix e in Eq. (^ at its experimental value. The other 
important quantity that needs to be calculated is the radiative width of the 0^ state, 
coming from the electric dipole (E2) decay into the 2^ state of ^^C This calculation 
involves the evaluation of the E2 operator between the initial 02" three-body scattering 
state and the final 2^ bound state [^. The proper three-body scattering-state treatment 



of the O2 initial state is not feasible in our approach for the time being, therefore we use 
a bound-state approximation to it. This is an excellent approximation for the calculation 
of r^ because the total width of the 02" state is very small (8.5 eV ll§]). The value of F^ 
is rather sensitive to the energy difference between the 02" and 2^ states, so we have to 
make sure that the experimental energy difference is correctly reproduced. 



2.3. Stellar-model calculations 

The composition of the interstellar material (ISM) is a mixture of ejecta from stars 
with different masses. At present it is not clear which type of stars contribute most of the 
^^C or ^^O to the ISM. Therefore, we performed stellar model calculations for a typical 
massive, intermediate-mass and low-mass star with masses 20, 5, and 1.3 Mq, respectively, 
including the calculated triple-alpha reaction rates. 

We used a modern stellar evolution code, which contains the most recent input physics 



27| . Up-to-date solar models can be calculated with this program [^ as well as the 
evolution of low-mass stars can be followed through the thermal-pulse phase of stars at the 
asymptotic giant branch ||29|. The nuclear network is designed preferentially to calculate 



the H- and He-burning phases in low-mass stars. Additionally, all basic reactions of C- 
and 0-burning are included, which may destroy the previously produced C and O in 
massive and intermediate-mass stars. 

Here, the stars are followed from the onset of H-burning until the third thermal pulse 
on the AGB, or until the core temperature reaches 10^ K in the case of the 20 Mq star 
(the nuclear network is not sufficient to go beyond this phase). 

Large portions of the initial mass of a star are returned to the ISM through stellar winds. 
Unfortunately, basically due to the simple convection model used in stellar modeling, the 
composition of the wind cannot yet be determined very accurately from stellar evolution 
theory. However, it is beyond the scope of the present investigations to determine how 
and when the material is returned to the ISM. Instead we examine how much C and O is 
produced altogether. 

For the 1.3 Mq star, which loses its envelope basically during the thermal-pulse phase, 
the maximum C and O abundances in the He-burning region have been extracted. Al- 
though the efficiency of the dredge-up of heavy elements to the surface is only badly 
known, it is independent of the nuclear physics, and hence should be similar in all models 
independent of the triple-alpha rate. By taking the maximum abundances in this region, 
we have a measure of how strongly the enrichment of the stellar envelope by C or O is 
altered by modifying the triple-alpha rate. 

In the 5 and 20 Mq stars further fusion reactions like C- and 0-burning take place. The 
dredge-up process of metal-enriched material under these circumstances is more compli- 
cated and more uncertain than in the 1.3 Mq star. The 20 Mq star finally even explodes 
in a supernova. Therefore, for the 5 and 20 Mq stars the total amount of C and O in the 
stellar interior is evaluated. 

3. VARIATIONS OF THE NUCLEON-NUCLEON FORCE 

In this section we slightly vary the strength and range parameters of the MN- and 
MHN-potentials and calculate the modified resonance energies e and gamma widths F^ 
of the 0^ state in ^^C With these values of e and F^ we recalculate the triple-alpha 
reaction rate r^.^ of Eq. (^. We find that the value of F-^ is very little changed by the 
small variations of the N-N interactions, leading to negligible changes in r^a- Thus, in 
the stellar model calculations we can fix F^ to its experimental value in all cases. The 
resonance energy e, however, is rather sensitive to variations in the N-N force, leading to 
large changes in the triple- alpha rate r^^. By making the N-N force weaker, the resonance 



energy e moves higher and the reaction rate rs^ decreases exponentially as can be seen 
from Eq. (^. The opposite is true for a stronger N-N force. 

We have done a series of tests in order to make sure that the main assumptions of 
our calculations remain valid while the interactions are varied. Here we mention only 
three main points. We checked that the use of Eq. (^ is justified for all the different 
resonance energies as well as temperatures occurring in our calculations. In all cases the 
triple-alpha reaction is dominated by the sequential process through the ®Be ground state 
and the 02" state of ^^C The amount of variations in the strength of the Coulomb and 
strong interaction are small enough to keep the experimental ^Be ground state and the 02" 
state of ^^C from becoming a broad state or a bound state. We also estimated the non- 
resonant contribution [^], without the 02" resonance in ^^C, to the triple-alpha rate and 
found that this contribution is about 7-13 orders of magnitude smaller than the resonant 
contribution in the stellar temperature range T k. (0.8 — 3) x 10^ K arising in our stellar- 
model calculations. Furthermore, the contribution of the next higher-energy 0^ resonance 
(at about 2.7 MeV above the three-alpha threshold |^2[) to the triple-alpha rate is smaller 
by about 40 orders of magnitude in the considered temperature range. Also, in the above 
temperature range the reaction will not proceed through the low-energy wing of the ^Be 
resonance 



30 



Some of the carbon, produced in the triple-alpha process, is further synthesized in the 
^^C(a,7)"'^^0 and ^^0(Q;,7)^°Ne reactions. The ^^0(a,7)^°Ne reaction is nonresonant, so 
variations in the strengths of the strong force can have only small effects on its reaction 
rate. The ^^C(a,7)^^0 process may look more dangerous because its cross section is 
strongly affected by subthreshold states in the ^^O nucleus p2[- However, if the N-N force 



is made weaker, then the subthreshold states become less bound, thereby enhancing the 
^^C(q;, 7)^^0 cross section. Therefore, in the case of a weaker force the small C/0 ratio is 
further decreased. An analogous reasoning holds for a stronger force. Thus, without doing 
any calculation for the ^^C(a, 7)^^© and ^^0(a, 7)^°Ne reactions with the modified forces, 
we can conclude that their effect would strengthen our hypothesis regarding carbon and 
oxygen production. 

There may be other windows for different values of the fundamental interactions in 
which an appreciable amount of carbon and oxygen could be produced. The window 
that seems to be closest to the values realized in our universe is the creation of oxygen 
and carbon in the big bang pT|. This possibility is prevented in our universe by the 



existence of the stability gaps at mass numbers A = 5 and A = 8. In order to make ^He 
or ^Li as well as ^Be become bound, we calculated that the strong forces would have to 
be increased by about 10%, a value that is about a factor 20 larger than the maximal 
variation investigated in this work. 

We multiply now in the MN- and MHN-potentials independently the strength param- 
eters Vofc as well as the range parameters tq^ for all repulsive and attractive terms V^k 
in Eq. (^) by the same factor p. This factor was set between 0.994 and 1.006 for the 
strength parameters and 0.997 and 1.003 for the range parameters of the MN and MHN 
potentials. With these values for the factor p we recalculate the resonance energies and 
the triple-alpha reaction rates. For each of these new reaction rates we then perform again 
the corresponding stellar model calculations. 

The resulting changes in the C and O abundances are shown (Fig. 0) with respect to 
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Figure 2. The change of the carbon (A) and oxygen (O) mass abundances {X) through 
variations of the strength (lower ordinate) and range parameters (upper ordinate) of the 
nucleon-nucleon potentials MN and MHN. The abundance changes are shown in panels a, 
b, and c for stars with masses of 20, 5, and 1.3 M©, respectively, in units of the standard 
values Xstand- The dashed curves are drawn to guide the eye. 



the case, where the standard value of the resonance energy has been used (i.e., with no 
variations of the strength or range parameters of the MN- and MHN-potentials). Because 
each shift in the resonance energy can be identified with a variation in the strength or 
range parameter of the corresponding N-N force, we scaled the lower and upper ordinate 
with variations in these quantities, respectively. 

We conclude that a change of more than about 0.5% in the strength parameters or 
0.2-0.3% in the range parameters of the MHN- or MN-potential would destroy either 
nearly all C or all O in every star. This implies that irrespective of stellar evolution the 
contribution of each star to the abundance of C or O in the ISM would be negligible. This 
corresponds to a fine tuning that is about two orders of magnitude better than obtained 



from the constraints of a bound deuteron [p^-|3^ or from the non-existence of a bound 
diproton or dineutron |135 |. 

In Fig. 1^ one can see that the size of the resulting change in the abundances with the 
variation of the range parameters r^k is about twice the size of the change coming from 
the variation of the strength parameters Vofc. This can be understood quite simply, be- 
cause we only considered small variations (~ 0.5 %) of the strength and range parameters, 
respectively. In this case the changes in the strength and range parameters in Eq. (B) 
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lead to a change in the spatial parts in the potential for 6rok <^ 1 
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where in the second line we used a mean potential range parameter r ^ rok in the first- 
order correction term. Therefore, in the above approximation a change in the range 
parameter rofc corresponds to a twice as large change, 6VQk ~ +26rok, in the strength 
parameter. 

4. ANTRHOPICALLY ALLOWED QUARK MASSES AND HIGGS PARA- 
METERS 

In this section we investigate the constraints on the anthropically allowed quark masses 



following largely the discussion already given in Ref. [36 



The one-boson exchange potential (OBEP) for the N-N force can be written as a sum 
of Yukawa terms 



VoBEp(r) 



-f 

J w 



2exp(— m^rj 



2 exp{-m^r) ^ exp^-m^r) 



(8) 



where the f 's and m's are the coupling constants and masses of the vr, a and u mesons, 
respectively. The OBEP-potential has a long-range tail from one-pion exchange, an at- 
tractive minimum due to the exchange of the hypothetical cr-meson and a repulsive core 
due to the cu-meson. In the following we will only consider the long-range part due to the 
pion, the one-pion exchange potential (OPEP), 
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For small changes of the pion mass, 
we obtain 
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where in the second line we used a mean potential range given by the pion mass r ^ l/m.,^ 
in the first-order correction term. 

The range defined by the inverse pion-mass is therefore also limited to a ±0.5% window. 
Here we assumed that the sensitivity of the carbon production to the potential strength 
is similar in the OPEP and MN/MHN cases, respectively. 

As long as the up- and down-quark masses are small compared to the QCD scale, the 



mass of the pion is well approximated by m.,^ oc J fTrijn^ + ^^d) [§3)01 • Therefore, the 
sum of the up- and down-quark masses rriu+d = m-a + m^ scales with the pion mass 



m^r as 0.55mu+d ~ (5m,r- Thus the anthropically allowed value of the sum of the up- 
and down-quark masses are fine tuned to approximately 1 %. Using the quark masses 
rriu = (4.88 ± 0.57) MeV and m^ = (9.81 ± 0.65) MeV the anthropically allowed sum of 
the up- and down-quark masses mu+d is constrained to within about 0.15 MeV. This is 



consistent with the result of Ref. [Q where a value in the order of 0.05 MeV is given. 

Constraints can also be obtained for the difference of the down- and up-quark masses, 
rrid-u = TTT'd — '^n- We do not consider possible variations of the electron mass as in 



Ref. |36|. The constraint for the upper limit of (5md_u follows from the restriction that 



the pp-fusion 

p-hp->d-Fe+-hz/ (12) 

is exothermic. This gives an allowed range of 5md-u < 0.42 MeV. The constraint for the 
lower limit of (5?7id-u follows from the stability of the hydrogen atom, i.e. that the proton 
capture of an electron 

p -|- e~ -^ n -|- z/ (13) 

is endothermic. This gives an allowed range of 5m^_^ > —0.782 MeV. The stability of 
the proton and deuteron give also constraints for the lower and upper limit of (5md-u- 
However, these constraints are less stringent than the ones obtained from the ones given 
above IBB . 



The resulting anthropically allowed strengths of the up- and down-quark masses are 
given by the central region of Fig. |^. As can be seen from this figure the anthropic fine 
tuning of the up- and down-quark masses is less than a few percent. The band around 
Sm^/m^ = —0.5Sm^/md determined by the production of carbon and oxygen is about 
four times narrower than the band around Smd/md = +0.56mn/rnu- That means that the 
constraints of carbon or oxygen production give a stronger fine tuning than the constraints 
determined from the reactions of Eqs. (O) and (|T^). If additional joint constraints that 



unification imposes on the fermion masses like the ratio md/rrie are taken into account [p7| , 
three-dimensional constraining plots for the up-quark, down-quark and electron masses 
can be derived. 



Similar considerations [^,^,^] can be carried out for the Higgs vacuum expectation 
value V ~ W — /i^/A oc /i and the Higgs mass parameter /x. Since the up- and down-quark 
masses scale simply with the Higgs vacuum expectation value, the Higgs vacuum expecta- 
tion value V and Higgs mass parameter /i are also fine tuned by anthropic considerations 
to approximately ±1 %. 

5. SUMMARY 

In this work we have shown that the strength or range of the strong force is fine tuned 
to ±0.5 % by the production of carbon or oxygen and therefore also carbon-based life. 
This leads to a fine tuning of the sum of the light quark masses as well as the Higgs 
vacuum expectation value and mass parameter to about ±1 %. 

One of the most fascinating aspects resulting from this is that life, stars, nuclei and 
elementary particles seem to be closely interwoven in our universe through this extreme 
fine tuning. It is not clear if the values of the fundamental parameters corresponding to 
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Figure 3. The central region shows the anthropically allowed values of variations of the 
up-quark mass 6mu/m^i and variations of the down-quark mass 6md/md. The follow- 
ing constraints are determined by the straight lines in the figure: (i) the pp-fusion in 
Eq. ([l^) is exothermic, (ii) the proton capture of an electron in Eq. ([I3| ) is endothermic, 
(iii) the carbon production is not suppressed, and (iv) the oxygen is production is not 
suppressed. These four constraints determine the upper and lower bounds of the bands 
around Smd/rrid = +0.56m^/m^^ and 6md/md = — 0.55mu/mu, respectively. 



this anthropic fine tuning will be calculable by a future Final Theory I^Q] or if some of 
the parameters of such a Final Theory will have to be chosen from a large or continuous 
ensemble [^. In any case it will be one of the challenges for a possible Final Theory to 
explain the anthropic fine tuning as described in this paper. 
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